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Tj- ; 1 Introduction 
00 

■ The investigation of Lie-point symmetry of differential equations is a well established and pow- 

CN ■ erful tool for constructing exact solutions of the equations under consideration [6l[T9 l [27 1 [29]. 

In fact, classical Lie-point symmetry group analysis is the only systematic method we know 
for deducing exact solutions (in particular, invariant, partially invariant and separation of vari- 
ables solutions) of nonlinear partial differential equations [SlITniESlI^SlEZlES]. However, for a 
lot of important applications differential equations because their classical Lie symmetry groups 
are rather trivial including at most space and time translations and scale transformations, the 
obtained solutions are not sufficient for our understanding the mechanics and mathematics of 
mathematical models. Thus, this stimulates the efforts devoted to generalization of Lie's original 
concept of symmetry in order to find more other types exact solutions. 

The first approach to these symmetries is conditional ones (called also nonclassical symmetries 
or Q-conditional symmetries), which was introduced by Bluman and Cole [5] in 1969. It consists 
in augmenting the original partial differential equations with invariant surface conditions, a 
system of first order differential equations satisfied by all functions invariant under a certain 
vector field. The latter is chosen as a classical infinitesimal symmetry for the augmented system. 
The basic equations for conditional symmetries are similar to Lie's determining equations except 
that they are nonlinear and less overdetermined. That is why one rarely succeeds in obtaining 
all possible solutions to the determining equations for conditional symmetries, especially in the 
case of second order equations. 

In contrast to classical Lie symmetry, another difference in the research of such two kinds of 
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symmetries is the procedure of deriving the determining equations. Namely, deriving systems 
of determining equations for conditional symmetries crucially depends on the interplay between 
the operators and the equations under consideration. For example, for the linear heat equation 
ut = Uxx, the determining equation are derived usually based on the condition that the general 
form of conditional symmetry operators Q = T{t, x, u)dt+(,{t, x, u)dx+r]{t, x, u)du, where (r, ^) / 
(0,0), were divided into two essentially different cases, i.e., the regular case r 7^ and the 
singular case r = 0. However, such partition of sets of conditional symmetries operators with 
the conditions of vanishing and nonvanishing coefficients of operators is not appropriate for all 
of differential equations. Therefore, investigation of algorithms for deriving the determining 
equations for condition symmetries are essential importance, which will provide some useful 
guidelines for simplifying the obtained determining equation and make its complete solving 
become possible. Recently, Popovych et.al give a detail investigation on this question and present 
a novel framework, namely singular reduction operators or singular reduction modules [8ll24). 
for finding an optimal way of obtaining the determining equation of conditional symmetries. As 
examples, properties of singular reduction operators of (l+l)-dimensional evolution equations 
and a specific wave equations have been studied p ^[30l432t [39 t HT] . However, for more general 
nonlinear wave equation, there exist no general results. 

In this paper we extend this new framework of singular reduction operators [21] to study a 
class of variable coefficient nonlinear wave equations which possesses three arbitrary functions 
/ = f{x),g = g{x) and h = h{x) and two power nonlinearities and of the form 



where f{x)g{x) 7^ 0, n and m are arbitrary constants, t is the time coordinate and x is the one- 
space coordinate. The linear case is excluded from consideration because it was well-investigated. 
We also assume the variable wave speed coefficient to be nonlinear, i.e. n 7^ 0. The case 
n = is quite singular and will be investigated separately. 

Many specific nonlinear wave models describing a wide variety of phenomena in Mechanics 
and Engineering such as the flow of one-dimensional gas, shallow water waves theory, longitudinal 
wave propagation on a moving threadline, dynamics of a finite nonlinear string, elastic-plastic 
materials and electromagnetic transmission line and so on, can be reduced to equation ([T]) (see [l] 
p. 50-52 and [2]). Classical Lie symmetries of various kinds of quasi-linear wave equations in two 
independent variables that intersect class ([I]) have been investigated in [3ll4l[6|[7tl9l[T0 l[T2Hl4l [T7l- 
[2 m [28 t[35]4371 H2]. Recently, we have present a complete local symmetry and conservation law 
classification of class ([TJ in |151I16]. Classical Lie symmetry reduction and invariant solutions 
of some variable coefficients wave models which are singled out from the classification results 
are also investigated. However, the nonclassical symmetries of class ([T|) remains open. In this 
paper, we will make a detailed investigation on this subjects so that we can use it to obtain some 
new non-Lie exact solutions. Below, following [24j we use the shorter and more natural term 
'reduction operators' instead of 'operators of conditional symmetry' or 'operators of nonclassical 
symmetry'. 

The rest of paper is organized as follows. In Section [21 known results on equivalence trans- 
formations and classical group analysis of the wave equations ([T|) are adduced in a form which 
is suitable for purposes of our investigations. Then in section [3l singular reduction operators, 
and in particular regular reduction operators classification for the class under consideration are 
investigated after some brief review of the notation of singular and regular reduction operator. 
Section H] contains nonclassical symmetry reduction of some nonlinear wave models including 
some truly 'variable coefficient' ones which are singled out from the class of (l+l)-dimensional 
wave equations. New non-Lie exact solutions of the models are constructed by means the re- 
duction. Conclusions and discussion are given in section [5j 
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2 Equivalence transformations and Lie symmetries 

The exhaustive result on classical group classification of class ([1]) is presented by the statements 
adduced below [15]. Using the direct method, we construct different kinds of equivalence groups 
including usual and generalized extended ones and discuss their structure (for the details about 
equivalence group can be seen in [21 ] [29 l [3"3 y 34 1 [38] ) . 

Theorem 1 . The usual equivalence group G~ for the class ([T]) consists of the transformations 



t = e2t + e3, x = X{x), u = e^u, 

f = eiele^^X-^f, g = eie^'^'^X^g, h = eie^'^X'^h, h = n, m 



where ej (j = 1, . . . ,4) are arbitrary constants, eie2e4 0, X is an arbitrary smooth function of 
X, X^ / 0. 

It is shown that class ([T]) admits other equivalence transformations which do not belong to 
G~ and form, together with usual equivalence transformations, an generalized extended equiva- 
lence group. Using the direct method [23], we can construct the following complete generalized 
extended equivalence group G~ of class ([T|). 

Theorem 2. The generalized extended equivalence group G~ of the class ([T]) is formed by the 
transformations 

t = e2t-\-es, x = X{x), u = Y{x)u, 

/ = eie2x-iy-"-2/, g = eoX,Y-^^-^g, h = eiX'^Y-^'-^-^h, h = n, m = m, 

where ej (j = 1,2,3) are arbitrary constants, ei€2 ^ 0, X is an arbitrary smooth function of x, 
Xx 7^ 0, the function Y{x) is determined by the formula 



1 



I (l-(n+l)F(x)) ^, f dx ^ 

[e^(^), n=-l J 9{x) 

Using the transformation 
~ f dx 

t = t, X = — -, u = u (2) 

J g{x) 

from theorem [H we can reduce equation ([TJ to f{x)uii = {u^Ux)x + h{x)u"^, where f{x) = 
g{x)f{x),g{x) = 1 and h{x) = g{x)h{x). (Similarly, any equation of form (1) can be reduced 
to the same form with f{x) = 1.) Thus, without loss of generality we can restrict ourselves to 
investigation of the equation 

f{x)utt = {u''ux)x + h{x)u^. (3) 

All results on symmetries, solutions and conservation laws of class ^ can be extended to class 
([1]) with transformations ([2]). 

Theorem 3. The usual equivalence group Gi for the class ([3]) with n ^ —1 is formed by the 
transformations 

u = [1 - (n + l)(e3X + £4)]"^^ = ((5Xa;)2^n =: Y{x)u, 
/ = efj^X^ h = 6 Xx h, n = n, m = m, 

where ej (j = 1, . . . ,4) and 6 are arbitrary constants, Jeiea / 0. 
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The above representation of the transformations from is very clumsily, but it shows that 
the family of transformations is continuously parameterized with the parameter n, including the 
value n = —1. In fact, we have lim„_j._i X{x) = f + limn_^_i Y{x) = e^3a;+e4_ Therefore, 

we can reformulate the above transformations as a more simple and clear formulae than in 
theorem [3l 

Theorem 4. The generalized equivalence group Gi for the class ([3]) with n / — 1 consists of 
the transformations 

i=eit + 62, X = f|2±f| =: u = e7X!^'u, 

3n+4 _ m + 3n+3 

/ = ef e?X^ /i = ef^+^+^X^ ^"+^ h, n = n, m = m, 

where ej {j = 1, . . . ,7) are arbitrary constants, eiey / and e^e^ — £465 = ±1. 
For n = —1 transformations from group takes the form 

i=eit + e2, x = e3X + e4, u = e5e'"^^u, 

/ = efeg^eg^e-^e^/, h = e^^e^'^e-"''^''h, h = n, m = m, 

where ej (j = 1, . . . ,6) are arbitrary constants, eie^e^ 7^ 0. 

Since the parameters n and m are invariants of all the above equivalence transformations, 
class ([1]) (or class ([3|)) can also be presented as the union of disjoint subclasses where each 
from the subclasses corresponds to fixed values of n and m. Thus the generalized equivalence 
group G"" of class ([T]) can be considered as a family of usual equivalence groups of the subclasses 
parameterized with n and m. Motivated by this, we can obtain wider equivalence groups for some 
of the subclasses apriori assuming the parameters n and m satisfy a condition. In particular, 
for the case m = n + 1 , we have 

Theorem 5. The class of equations 

f{x)utt = {g{x)u^u^)^ + /i(x)n"+i (4) 
admits the equivalence transformation group G~^„_j_| consisting of the transformations: 

i=e2t + e3, i = X(x), u = Y{x)u, / = eie^X-iy-^-^, 

g = eiX,Y-^^-^g, h = ei[h - Y^+HY-^''+^^YMX-^y'^'''% (n = n), 

where X andY are arbitrary functions ofx, ej (j = 1,2,3) are arbitrary constants, €i€2XxY 7^ 0. 

If we further gauge g with the condition g = 1 and impose the restrictions eiX,^. = y2n,+2 
parameters of G^^^^i, we can obtain the following equivalence group for the class of 

equations 

f{x)uu = {u''ux)x + h{x)u^+\ (5) 
Corollary 1. The equivalence group Gl,^^^_^_l of class ([5]) is formed by the transformations: 
i=e2t + e3, x = X{x), u = Y{x)u, 

/ = e2g2y-3n-4j^ /i = ef [/i - Y^+i (y-("+2)y^,)^]y-4n-4;j^ (n = n), 

where X and Y are arbitrary functions satisfying the condition eiX^ = y2n+2^ ^. _ 1^2,3) 
are arbitrary constants, eie2Y 7^ 0. 
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Using the above-mentioned results, we obtained the classical symmetry classification for the 
class ([3]) in [15] . 

In table dl we list tuples of parameter-functions /(x) and h{x), the constant parameters m 
and n and bases of the invariance algebras in all possible inequivalent cases of Lie symmetry 
extension. The operators from table [1] form bases of the maximal Lie invariance algebras iff 
the corresponding values of the parameters are inequivalent to ones with more abundant Lie 
invariance algebras. In cases [2] and [3] we do not try to use equivalence transformations as much 
as possible since otherwise a number of similar simplified cases would be derived, see footnote 
after the table. 



Table 1. Results of group classification of class ((3]) 



N 


n 




h[x) 


Basis of A"""" 


General case 


f 


V 




V 


V 


dt, 


2 


V 






h'(x) 


dt, \{d+2h- pn)tdt + [a{n+ l)x^ + hx + c]da: 
+ {ax + p)udu, 


3 


V 




1 


e 


dt, dx, (1 — m)tdt + (1 + n — m)xdx + 2udu 










m — 


l,h^O, {h/f)x=0 


4 






fix) 


ef 


dt, n[a{n + l)x'^ + bx + c]dx + {nax + 2b + d)udu 


5 




4 
3 


1 


e 


dt, dx, nxdt + 2udu 


6 


-4 




/'W|n=-4 


ef 


dt, e^\dt + udu), e-^\dt ~ ua„),4(-3ax2 + 62; + c)dx 
+ {'lax - d - 2h)udu 


7 


4 
3 




1 


e 


dt, dx, 2xdx — 3udu,x^dx — Sxudu 


8 


-4 




f 


1 


dt, dx, e^^dt + udu), 2xdx - it9„, e"2'(9t - udu) 


9 


-4 




f 


-1 


dt, dx, cos{2t)dt - sm{2t)udu, 2xdx - udu, 
sm{2t)dt + cos{2t)udu 


m = n + lor/i = 


fO 


/-4 




V 


V 


dt , ntdt — 2udu 


ff 




4 
3 


1 




dt, tdt + xdx, ntdt — 2udu, 


12 




4 
3 


f 


e 


dt, dx, ntdt - 2udu 


13 




4 
3 


1 





dt, dx, tdt + xdx, ntdt — 2udu 


U 


4 
3 






Q 


dt, tdt + 2dx, 2tdt+Sudu 


15 


4 
3 




1 





dt, dx, tdt + xdx, 2tdt + Sudu, x^dx — Sxudu 


16 


-4 




f 





dt, dx, tdt + xdx, 2tdt + udu,t^dx + tudu 



Here a is arbitrary constant, e = ±1, 

and it can be assumed up to equivalence with respect to G\ that the parameter tuple (a, b, c, d) takes only the 
following non-equivalent values: 

(£,0,1,0) if n=-l or (1,0, if n ^ -1, (0, l,0,d')> (0,0,1,1), 

where d' is arbitrary constant. In all the cases we put g(x) = 1. In case [10] the parameter-functions / and h 
can be additionally gauged with equivalence transformations from Gi^m=n+i- For example, we can put / = 1 if 
n 7^ —4, —4/3 and / = otherwise. 

Remark 1. It should be emphasized that we adduce only the cases of extensions of maximal 
Lie invariance algebra, which are inequivalent with respect to Gi \i m ^ n + \ and with respect 
to G-^^ 

=n+i if ?Ti — n -|- 1 or /i — in table [H What's more, right choice of a gauge of the 
parameter tuple (/, 5, h) is another crucial point of our investigations. In fact, the major choice 
have been made from the very outset of classification when the parameter-function g was put 
equal to 1. It is the gauge that leads to maximal simplification of both the whole solving and 
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the final results. Although all gauges are equivalent from an abstract point of view, only the 
gauge g = 1 allows one to exhaustively solve the problem of group classification of class ([1]) 
with reasonable quantity of calculations. Even after the other simplest gauge / = 1 chosen, 
calculations become too cumbersome and sophisticated. 



3 Nonclassical symmetries 

In this section, we first review some necessary definitions and statements on nonclassical symme- 
tries [TTl|25ll3Qlll3] and singular reduction operator [Mj. Then perform a detailed investigation 
on the reduction operators of class ([3]). 

Consider an rth order differential equation C of the form L{t,x,U(^j.)) — the unknown 
function u of the two independent variables t and x, where L = L[u] = L{t,x,u^j.^) is a fixed 
differential function of order r and it(^) denotes the set of all the derivatives of the function u 
with respect to t and x of order not greater than r, including u as the derivative of order zero. 

In order to discuss the conditional symmetries of equation C, we will first treat equation C 
from a geometric point of view as an algebraic equation in the jet space J*" of order r and is 
identified with the manifold of its solutions in [27] : 

£ = {(t,X,U(^r)) £ J^\L{t,X,U(^r)) = 0}- 

Let Q denote the set of vector fields of the general form 

Q = T{t,x,u)dt + ^{t,x,u)d^ + ri{t,x,u)du, (r,^)/(0,0), (6) 

which is a first-order differential operator on the space x with coordinates t,x, and u. 
Then all functions invariant under Q and only such functions satisfy a first order differential 
equation 

Q[u] := TUt + (u:c - r] = (7) 

called the the characteristic equation (also known as invariant surface condition). 

Denote the manifold defined by the set of all the differential consequences of the characteristic 
equation Q[u] =0 in J''' by Q{r)) i-e., 

Q(^) = {(t,x,n(^)) G r\DfD^Q[u] = 0, a,/3 G NU{0}, a + < r}, 

where Dt = dt + Ua+i^pdu^ p and Dx = dx + Ua^p+idu^ ^ are the operators of total differentiation 
with respect to the variables t and x, and the variable Ua^p of the jet space J*" corresponds to 
the derivative 

Denote also by Q(^y) the standard rth prolongation of Q to the space J^: 

0<a+/3<r 

Definition 1. The differential equation C is called conditionally invariant with respect to the 
operator Q if the relation 



Q(.r)[L(t,x,n(r))] 



= (8) 



holds, which is called the conditional invariance criterion. Then Q is called conditional symmetry 
(or nonclassical symmetry, Q-conditional symmetries or reduction operator) of the equation C. 
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We denote the set of reduction operators of the equation L by Q(£) which is a subset of 
Q. Any Lie symmetry operator of L belongs to Q(/2). Sometimes Q(>C) is exhausted by the 
operators equivalent to Lie symmetry ones in the sense of the following definition. 

Definition 2. Two differential operators Q and Q in Q are called equivalent (Q ^ Q) if 
they differ by a multiplier which is a non- vanishing function of t,x and u . Q = XQ, where 
A = X{t,x,u),X ^ 0. 

Consider a vector field Q in the form ([6]) and a differential function L = L[u] of order ord L = r 
(i.e., a smooth function of variables t,x,u and derivatives of u of orders up to r). 

Definition 3. The vector field Q is called singular for the differential function L if there exists 
a differential function L = L[u] of an order less than r such that L\q^^^ = L\q^^^ . Otherwise Q 
is called a regular vector field for the differential function L. If the minimal order of differential 
functions whose restrictions on Q^^) coincide with L\q^^^ equals k (k < r) then the vector field 
Q is said to be of singularity co-order k for the differential function L. The vector field Q is 
called ultra- singular for the differential function L if L\Q(r) = ^■ 

Definition 4. A vector field Q is called weakly singular for the differential equation C : L[u] = 
if there exists a differential function L = L[u] of an order less than r and a nonvanishing 
differential function A = X[u] of an order not greater than r such that L\q^^^ = XL\q^^^ . 
Otherwise Q is called a weakly regular vector field for the differential equation C. If the minimal 
order of differential functions whose restrictions on Q(^) coincide, up to nonvanishing functional 
multipliers, with L\qi^^^ is equal to k{k < r) then the vector field Q is said to be weakly singular 
of CO- order k for the differential equation C. 

Definition 5. A vector field Q is called a singular reduction operator of a differential equation 
£ if Q is both a reduction operator of £ and a weakly singular vector field of £. 

After the factorization of the reduction operator under the usual equivalence relation of re- 
duction operators in Definition [2] to singular and regular cases, the classification of reduction 
operators can be considerably enhanced and simplified by considering Lie symmetry and equiv- 
alence transformations of (classes of) equations. 

Lemma 1. Any point transformation of t, x and u induces a one-to-one mapping of Q into 
itself. Namely, the transformation g :i = T(t,x,u), x = X{t, x,u),u = U{t,x,u) generates the 
mapping g* : Q ^ Q such that the operator Q is mapped to the operator g^Q = fd^ + ^dx + fjdu, 
where f{i,x,u) = QT(t,x,u), S,{i,x,u) = QX{t,x,u), fj{i,x,u) = QU{t,x,u). If Q' ^ Q then 
g*Q' ~ g*Q- Therefore, the corresponding factorized mapping gf ■ Qf ^ Qf also is well defined 
and bijective. 

Lemma [T] results in appearing equivalence relation between operators, which differs from 
usual one described in Definition [2j 

Definition 6. Two differential operators Q and Q in Q are called equivalent with respect to a 
group G of point transformations if there exists g G G for which the operators Q and g^:Q are 
equivalent. We denote this equivalence by Q ~ Q mod G. 

The problem of finding reduction operators is more complicated than the similar problem 
for Lie symmetries because the first problem is reduced to the integration of an overdetermined 
system of nonlinear PDEs, whereas in the case of Lie symmetries one deals with a more overde- 
termined system of linear PDEs. The question occurs: could we use equivalence and gauging 
transformations in investigation of reduction operators as we do for finding Lie symmetries? 
The following statements give the positive answer. 
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Lemma 2. Given any point transformation g of an equation C to an equation C, g^, maps Q{C) 
to Q{C) bijectively. The same is true for the factorized mapping gj from Qf{C) to Qf(C). 

Corollary 2. Let G he the point symmetry group of an equation C Then the equivalence of 
operators with respect to the group G generates equivalence relations in Q{C) and in Qf{C). 

Consider the class £1^ of equations Cg: L(t, x, U(j.^,6) = parameterized with the parameter- 
functions 9 = 6{t, X, U(^r))- Here L is a fixed function of t, x, U(^r) and 9. The symbol 9 denotes the 
tuple of arbitrary (parametric) differential functions 9{t, x, = {9^{t, x, 9^{t, x, 
running through the set S of solutions of the system S{t, x, 9(^g^{t, x, it(r))) = 0. This system 
consists of differential equations on 9, where t, x and play the role of independent variables 
and 9^q) stands for the set of all the derivatives of 6 of order not greater than q. In what follows 
we call the functions 9 arbitrary elements. Denote the point transformation group preserving 
the form of the equations from £1^ by 

Let P denote the set of the pairs consisting of an equation Cg from and an operator Q 
from Q(£g). In view of Lemma [2l the action of transformations from the equivalence group G~ 
on C\s and {Q{C9)\9 € S} together with the pure equivalence relation of differential operators 
naturally generates an equivalence relation on P. 

Definition 7. Let 9,9' e S, Q e Q{Ce),Q' G Q[Ce'). The pairs {Ce,Q) and {Ce',Q') are 
called G~-equivalent if there exists a g £ G~ such that g transforms the equation Cq to the 
equation£e', and Q' ^ g^Q. 

The classification of reduction operators with respect to G~ will be understood as the clas- 
sification in P with respect to this equivalence relation, a problem which can be investigated 
similar to the usual group classification in classes of differential equations. Namely, we construct 
firstly the reduction operators that are defined for all values of 9. Then we classify, with respect 
to the values of 9 for which the equation Cq admits additional reduction operators. 

In what follows we use above-mentioned method to investigate reduction operators of the 
(l-l-l)-dimensional variable coefficient nonlinear wave equations For convenient, we rewrite 
it as the form 

f{x)uu - Kn,.)x - Hx)u"' = 0. (9) 
3.1 Singular reduction operators 

Using the procedure given by Popovych et al. in |24j, we can obtain the following assertion. 

Proposition 1. A vector field Q = T{t,x,u)dt + ^{t, x,u)dx + r]{t, x,u)du is singular for the 
differential function L = f{x)utt — {u^Ux)x — h{x)u"^ if and only if ^^/(x) = t'^u^ . 

Proof. Suppose that r 7^ 0. According to the characteristic equation rut + (,Ux — rj = 0, we can 
get 

r] C 

ut = Ux, 

T T 

Utt = {-)t - [-)tUx + [{-)u - (.-)uUx\[ Ux) - (,-)[[-)x + {-)uUx - {-)xUx 

T T T T T T T T T T 



XX I ■ 
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Substituting the formulas of Uu from above formulaes into L, we obtain a differential function 
L = [f{x){-f - u^]u,, + f{x){{\ - {^)tu., + [C^)u - {-)uuM- - -n.) 



-(-)[(-). + {-)uu. - (-)xUx - {-)uul] \ - nu^-^ul - hu^u 



T T T T T 



X ■ 



According to the definition [3] of singular vector field, we have ord L < 2 if and only if f{x){-)'^ — 
-u" = 0. ^ □ 

Therefore, for any /, h, n and n with /u" > the differential function L = f(x)utt — {u^Ux)x — 
h{x)u^ possesses exactly two set of singular vector fields in the reduced form, namely, S = 
{dt + \/u^ / fdx + fidu} and S* = {dt — \JuP- 1 jdx + where t) = ^. Any singular vector 

field of L is equivalent to one of the above fields. The singular sets are mapped to each other 
by alternating the sign of x and hence one of them can be excluded from the consideration. 

Proposition 2. For any variable coefficient nonlinear wave equations in the form ([9]) the dif- 
ferential function L = f{x)utt — {vP'Ux)x — h{x)u"^ possesses exactly one set of singular vector 
fields in the reduced form, namely, S = {dt + \/u"/ fdx + fjdu}- 

Thus taking into accountant the conditional invariance criterion for an equation from class 
([9]) and the operator dt + \/u'^/ fdx + r]du , we can get 

Theorem 6. Every singular reduction operator of an equation from class ([9j) is equivalent to 
an operator of the form 



Q = dt + \/u'^/ f{x)dx + r]{t, x, u)du, 

where the real-valued function r]{t,x,u) satisfies the determining equations 

{-l/2nhu^-^ - 2fr]tu - 2/r?r?„,)7^ + (3/4/2// _ i/2/,,)(nV/)3/2 
+ (-l/4n2r/2ti"-2 _ n7?r?„u"-i + l/2n7?2u"-2)(n«//)-i/2 _ nry^^'^-i 

-l/2nrju^-^fx/f - 2r/,„n" = 0, (10) 
-r]xxu"- - nr]r]xu'^~^{u'^ / f)-^/"^ + h-quu"" + 2friritu + frfr]^ 



luu 



-mhvu"'-^ + frjtt + {hfx/f - hx)u'^^/^ = 0. 
3.2 Regular reduction operators 

The above investigation of singular reduction operators of nonlinear wave equation of the form 
([9]) shows that for these equations the regular case of the natural partition of the corresponding 
sets of reduction operators is singled out by the conditions ^ ^ -^^Ju^ j fr . After factorization 
with respect to the equivalence relation of vector fields, we obtain the defining conditions of 
regular subset of reduction operator: r = 1,^ 7^ zby^n"-//. Hence we have 

Proposition 3. For any variable coefficient nonlinear wave equations in the form ([9]) the dif- 
ferential function L = f(x)utt — {u^Ux)x — h{x)u'^ possesses exactly one set of regular vector 
fields in the reduced form, namely, S = {dt + ^dx + fjdu} with ^ / ±y^u^/ f. 

Consider the conditional invariance criterion for an equation from class ([9|) and the operator 
dt + £,it,x,u)dx + r]{t,x,u)du with ^(t, x,u) / it y^u"//, we can get the following determining 
equations for the coefficients ^ and rj: 

e„ = 0, 2f^t - nrju^-' + {2^x + Ux/ f)u^ = 0, 

{2n£,x - nrju + nS,fx/ f)vJ'~^ + (nr] - n^r])u'''-^ - rjuuu'^ + /?^7?«« = 0, . . 

2fCt^x - 2fCtVu - 2nr]xu^-^ - f^tt - ^fii^iiuu - ViVtu + {^xx - 2vxu)u^ = 0, ^ ^ 

i^hfx/f - ^hx + hr]u)vr + frfrjuu + 2fi]i]tu - 2f^tVx - VxxU"^ + fvtt - mhr]u"^-^ = 0. 
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From the first two equations of system (jlip . we have 

n n J 

Substituting these expression into the last three equations of system we have the fohowing 
assertion. 

Theorem 7. Every regular reduction operator of an equation from class ([9]) is equivalent to an 
operator of the form 

Q = dt + Cit,x)d.^+r]{t,x,u)du with 7?(t,x,n) = -2/^4^1-" + -(2ex+?^)n, (12) 

n n J 

where the real-valued function ^(t, x) satisfies the overdetermined system of partial differential 
equations 

2(l-n)/6 = 0, 2(l-n)mi = 0, 8(1 - n)/^^^^ ^ 0, 
4(1 - n)[(/6)2 - f^U^^. + i^) + f^itt] = 0, 

4,-2(i + i)(2e. + e^), = 0, 

Viti. - lfU2i:, + e^) - 4(/6). - fitt - Ifmx + i^)t - ^(1 - n)(/6)x = 0, 
{ihif - ih^ + i/i(l - m)(2C, + iff))u"^ + 1(1 _ n - m)fhitu^-^ (13) 
-^{l - n)f\itfu^-^- + - n)[f\t{fit)t - fHm. + ^^)]ni-2" 
+ [lf{fit)tt - ^fUm. + ^f%{2^. + i^)t - ^1 - n)f%{2^^ + iiff 

+^(1 - n)f{fit)t{2i, + iff)y-^ + [^/(2e. + i^fm., + ^^)t 
+^/(2ex + ^^)tt - ^/6(2ex + e^)x - Uf^th.]u - ^(2^. + ^fun^+' = o. 

Solving the above system with respect to the coefficient functions S,, f and h under the 
equivalence group Gi, we can get a classification of regular reduction operator for the class ([9]). 
It is easy to know that some of the regular reduction operator are equivalent to Lie symmetry 
operators, while some of are nontrivial. Below, we give a detailed investigations for these cases. 

In fact, the first three equations of system (fT3|) implies there are two cases should be consid- 
ered: n 7^ 1 or not. (It should be noted that ^ = should be exclude from the consideration 
because it leads to r/ = 0). 

Case 1: n ^ 1. In this case, we have = 0. Thus system ([T3|) can be reduced to 
(3n + 4)C,, + 2(n+l)(e^), = 0, 

(e/i^ - ^h., + i(l - m)h{2^^ + e^))tx- - i(2e.. + e^)....n"+i = 0. ^ ^ 

Thus, there are two cases should be considered: m 7^ n + 1 or not. 

Case 1.1: For m 7^ re + 1, from the second equation of ()14p we obtain 

ih^-f -ih^ + -{l- m)h{2i^ + ify) = 0, (2e. + C^)^^ = 0. (15) 
/ n f ^ f '''' 

The first equation of suggests that (3re + 4)^^: + 2(n + l)^y is independent of the variable 

X, so there exists a constant r such that (3re + 4)^x + 2(n + 1)^^ = nr. The second equation of 

(jlSp suggests that there exist two constants a and b such that 2^^ ~ ^^"^ + '^^- By solving 

the last two equations we obtain 

f 

= 2(re + l){ax + b)-r, = 2r - (3n + A){ax + b), 
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which together with the first equation of (jl5p imply 

^ = a(n + + [2b{n + 1) - r]x + s, 

/■ 2r-(3n + 4)(ax + b) . 

, /■ 2r - (m + 3n + 3)(ax + 6) . 
(xj-exp(y _^ ;L)2;2 + [26(n + 1) - r]x + s 

where a, 6, r, s are arbitrary constants. Thus, the corresponding regular reduction operator has 
the form 

Q = dt + [a{n + l)^^ + (26(n + 1) - r)x + s]d^ + {ax + b)udu, 
which is equivalent to Lie symmetry operator. 

Case 1.2: m = n -\- 1. In this case, system (jl4p can be rewritten as 

La 

J'-^ (17) 



(3n + 4)C^^ + 2(n + l)(e^)^ =0, 



Integrating these two equations with respect to functions f{x) and g{x), we can obtain 



371+4 . /■ 1 . £ — 2££^.^. - 

/(x) = |er^exp(ry-dx), /^(x) = ^^^^^1^ 



where p, r are arbitrary constants, £ is an arbitrary smooth function and n ^ —1. In addition, 
r] = ;^(2^x + = + 2ra+2 )^' Thus, we have a nontrivial regular reduction operator 

g = + i{x)d^ + [(^ + ^-ii-^)n]a„, n / -i. (is) 

It should be noted that for n = — 1 the reduction operator is also equivalent to Lie symmetry 
operator. 

Case 2: n = 1. In this case, we have = 2/£t + (2^^^ + £^)u. Thus, system can be 
reduced to 

(TCx + 4^=^)^ = 0, 2[£^ + 2(e + l)^]6 + 4(e + l)6x+6i = 0, 
[2/(2C. + e^)(2e, + £^)t + /(2e.. + iif)tt - 2fCt{2^^ + C^)x - 2(/6)xx] u 
-(2ex + ^^)xxu^ + 2f{2ititx + 6«) + [e/i^ - e/ix + (1 - m)h{2i^ + e^)]n" 
-2mhfCtu"'~^ = 0. 

After some brief analysis, we find that there are five cases should be considered. 
Case 2.1: m = 0. In this case, the third equation of (|19p implies 



(19) 



2h^ff - h,^ + 2h^, + if^t^tx + 2f%u = 0, 



2/(2^. + ^fmx + cf)t + f{2ix + if)tt - 2fit{2ix + if)x - 2{fit)xx = o, (20) 

(2^3:: + (,~^)xx — 0. 

From the last equation of system (|20p we can know that there exist two functions a{t) and b{t) 
such that 2^x + ^.y- = a{t)x + b{t). On the other hand, the first equation of (|19p implies there 

exists a function c{t) such that 7^a; + 4^^ = c(t). Solving the last two equations gives 

f 

= 4a{t)x + 4b{t) - c{t), £y = -7a{t)x - 7b{t) + 2c{t), 
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from which we can get 



^ = 2a{t)x^ + 4:b{t)x - c{t)x + d{t), 

f(x)-e.J f^J^(th^m±M^i^] (21) 
/ [X) -exp^j ^^^^^^2 + 45(^)2; _ + d{t) ^ ■ 

where d{t) is an arbitrary functions. Since y is independent of t, we see that 
-la{t)x - 7b{t) + 2c{t) 



-2a{t)x^ + 4b{t)x - c{t)x + d{t) 
which leads to 



0, 



'U[a{t)b'{t) - a'{t)b{t)] + 3[a' {t)c{t) - a{t)c' (t)] =0, 

< [b{t)c'{t) - b'{t)c{t)] + 7[a{t)d'{t) - a'(t)d(t)] = 0, (22) 
2[c'{t)d{t) - c{t)d'{t)] + 7[b{t)d'{t) - b'{t)d{t)] = 0. 

Now, we multiply both sides of the second equation of ([T9|) by and substitute (f2T|) into it, 
then simplify the equation and compare the coefficient of x*(i = 0, 1, . . . , 5) to obtain 

'a'(t) = 0, 

a^[-4b'{t) + c'{t)] = 0, 

a[8c{t)c'{t) + 20ad'{t) + c"{t) - 4b"{t) - ?>2c{t)b'{t) + ll2b{t)b'{t) - 2Sb{t)c' {t)] = 0, 
-Ac"{t)b{t) + 32ac{t)d'{t) - 80ab'{t) + 20ac'{t) - 4b"{t)c{t) - 120a6(i)d'(i) 
+c"{t)c{t) + 2d"{t)a + 2c^{t)c'{t) - 12b{t)c{ty{t) + m{tfc'{t) - 8c\t)b'{t) 
+m{t)c{t)b' (t) + 16b"{t)b{t) - 64b'^{t)b'{t) - 16ad{t)b'{t) + 4ad{t)c'{t) = 0, (23) 
[-20ad'{t) + m{t)b'{t) + 2c{ty{t) - I2b{t)c'{t) - 8c{t)b'{t) + Ab"{t) - c"{t)]d{t) 
-80b'^{t)d'{t) + 4:d"{t)b{t) - d"{t)c{t) - 4c(t)c'(t) + 446(t)c(t)d'(t) - 28ad'{t) 
-m{t)b'{t) - Qc^{t)d'{t) + 12b{t)c'{t) + I6c{t)b'{t) = 0, 

[166'(t) - Ac'{t)]d^{t) + [-20b{t)d'{t) - 4c'(t) + 6c{t)d'{t) + d"{t) + lQb'{t)]d{t) 
+8c{t)d'{t) - 28b{t)d'{t) = 0. 

Note that ^ is assumed not to be identical with zero, after some simple but lengthy compu- 
tations, we find that systems ([2^ and ([25]) can be reduced to: 

a'(i) = 0, &'(t)=0, c'(i) = 0, (i'(t) = (24) 



or 



or 



or 



a = 0, c(t) = 46(t), (i(t) = g6(t), g6"(t) + Aqb{t)b'{t) + 4&'(t) = (25) 



a = 0, 2c(t) = 76(t), 6"(t) = -36(t)6'(t), b{t)d'{t) + 2b' {t){d{t) + 1) + d"{t) = (26) 



a = 0, c{t) = 3b{t), d = qb{t), b" (t) + 2b{t)b' (t) = 0, (27) 
where q is an arbitrary constant. 
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Case 2.1a: If system ()24p is satisfied, then = 0, the second equation of (j20p is an identity. 
The expression (j2ip can be rewritten as 



„, , , , 7ax-n + 2c ^ \ (28) 
/ (x) = exp / - — 7p——r — d X , 



^ = 2ax^ + Ahx — cx + d. 

—lax 

2ax^ + Ahx — cx + d 

where a, b, c and d are arbitrary constants. The first equation of ()20p is reduced to 

Substitute the expression of ^ and /(x) into it and integrate both sides to obtain 

, , . ( f -Qax - 66 + 2c , \ 
/i X = exp / ^ ; :dx . 

^ ^ 2ax2 + 462; - cx + d / 

In addition, ry = 2/^t + (2^a; + i^)u = {ax + 6)n. Therefore, we have 

^ = 2ax^ + Ahx — cx + d, 
7] = {ax + 6)n, 

/(x)=expf / -7"^-'^ + ^^ dx), (29) 



2ax^ + 46x — cx + d 
, , , / /■ -6ax - 66 + 2c . \ 

^^^)=^"Ply 2ax2 + 46x-cx + /"> 

where a, 6, c, d are arbitrary constants. Thus, the corresponding regular reduction operator has 
the form 

Q = dt + (2ax^ + 46x — cx + d)dx + {ax + b)udu, 
which is equivalent to Lie symmetry operator. 



Case 2.1b: If system (j25|) is satisfied, then the expression ([21]) can be rewritten as 

^ = qb{t), /(x)=exp(^). (30) 

Hence, = 0, k = b{t). Substituting these formulaes into the second equation of (pOj) we obtain 

qb"{t) + 2qb{t)b'{t)-2b'{t) =0. 

Combine it with the fourth equation of ([25]) to get b'{t) = 0. Hence a{t) , b{t) , c{t) , d{t) satisfy 
system (|24|) . and the solution is included in the case 2.1a. 



Case 2.1c: If system (I26p is satisfied, then the expression (I2ip can be rewritten as 

C = ^b{t)x + d{t), /(x) = l modG^. (31) 

Substitute it into the second equation of ()20p to obtain 2b{t)b'{t) + 6"(t) = 0. Combine it 
with the third equation of (f26l) to get 6'(t) = 0. Substitute it into the fourth equation of (f26|) 
to obtain bd'{t) + d"{t) = 0, which implies d{t) = 716"''* + 70, where 71 and 70 are arbitrary 
constants. Therefore ^ = |x + 7ie~^* + 7o- Substitute it into the first equation of (pUt) to obtain 

2-fi{hx + 26^)e"^* + bxhx + 2johx - 2bh = 0. 
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Since /i, 71,70 are independent of t, the preceding equation suggests that 

27i(/i^ + 263) = 0, bxh,^ + - 2bh = 0, 

which leads to 6 = 0. Therefore, we have 

(^ = dit + dQ, T] = 2do, / = 1, h = ho mod 
where di,do, /iq are constants. Thus, the corresponding regular reduction operator has the form 

Q = dt + {dit + do)dx + 2diudu, 
which is equivalent to Lie symmetry operator. 



Case 2. Id: If system (|27|) is satisfied, then the expression (|2T|) can be rewritten as 

^ = b{t){x + q), f{x) = ^— modG^. (32) 

X + q 

Substitute it into the second equation of (|20p to obtain 2b{t)b' {t) + b" (t) = 0, which is equivalent 
to the fourth equation of ([27]), and which leads to 2b''^{t) + b"'{t) = -2b{t)b"{t). Substitute 
(|32p into the first equation of ()20p to obtain 



K = 2f[2b'\t) + b"\t)]/b{t) = 2f[-2b{t)b"{t)]/b{t) = -4fb"it). 

Since h and / are independent of t, there is a constant r such that b"{t) = r. It follows that 
there exist constants s and w such that b{t) = rt^ /2 + st + w. Substitute it into the fourth 
equation of ([27|) to obtain 

rH^ + 2,rst^ + 2{wr + s^)i + 2ws + r = 0. 

Then r = and s = 0. Hence b{t) = w, = 0, the solution is included in the case 2.1a. 



Case 2.2: When m = 1, the third equation of (jl9p implies 
.Jk^eh _i_ fi'ot _i_ f ikv. _L o f rot J- cJk\(oe j^eL 



-2/6(2^, + e^), - 2ifCtU = 0, .33^ 



Similar to the case of in = 0, from the third equation of (jSSp and the first two equations of 
(fT9l) we get the expression of ^ and f{x) as stated in (i2T]) . where a{t) , b{t) , c{t) , d{t) satisfy the 
condition (l2ll or (1251) or (|26D or (127 



Case S.Sa.' If system is satisfied, then = 0, the second equation of (p3|) is an identity. 
The expression (i2T]l can be rewritten as (f28]l . The first equation of (p3]) is reduced to hx/h = 
fx/ f, which leads to h{x) = ef{x) (e = ±1) mod . In addition, ry = 2/^t + {2^x + i^)u = 
[ax + b)u. Thus, we have 

^ = 2ax^ + 4hx - cx + d, 
7] = (ax + b)u, 

= ( / 2ax2T4J-tJ+d ^ ^) ' ^''^ 
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where a,b,c,d are arbitrary constants and e = ±1. Thus, the corresponding regular reduction 
operator has the form 

Q = dt + {2ax^ + 46a; — cx + d)dx + {ax + b)udu, 

which is equivalent to Lie symmetry operator. 

Case 2.2b: If system (I25p is satisfied, then the expression (I2ip can be rewritten as ()30p . Hence, 
Cx = 0, k = b{t). If b'{t) = 0, then a{t) , b{t) , c{t) , d{t) satisfy system (p^ . and the solution is 
included in the case 2.2a. We suppose that b'{t) / 0. From the second equation of (f33l) we see 
that h = fCttt/Ct- Substitute it into the first equation of ([55]) to get fktt + 2fkkt — 2fxx(,t = 0- 
Further it can be reduced to qb"{t) + 2qb{t)b'{t) — 2b'{t) = 0. Combine it with the fourth 
equation of ([251) to get b{t) = —3/q which is contradict to the hypothesis b'{t) / 0. 



Case 2.2c: If system (|26|) is satisfied, then the expression (|21|) can be rewritten as (|31|) . If 
b'(t) = d'{t) = 0, then a{t) , b{t) , c{t) , d{t) satisfy both systems (p^ and ([26]) . and the solution 
is included in the case 2.2a. We suppose that b'^{t) + / 0. Substitute §Bj into the first 

equation of ([33|) to obtain 

b{t)x + 2d{t) ^ ^ 

Substitute (f3T|) into the second equation of (f33l) to obtain 

Then 

_ 2[b"'{t)d'{t) -b'{t)d"'{t)] 
~ [6'(t)x + 2d'(t)]2 ■ 

Substituting it into ([35]) yields 

b"'{t)d'{t) -b'{t)d"'{t) _ b"{t) + 2b{t)b'{t) 
[b'{t)x + 2d'{t)\'^ ~ b{t)x + 2d{t) 

Compare the coefficient of to obtain b''^{t)[b" [t) + 2b{t)b'{t)\ = 0. Substitute the third 
equation of (f26]l into it to obtain b{t)b'^{t) = 0, hence b'{t) = 0. Thus the fourth equation of 
(j26p can be reduced to bd'{t) + d"{t) = 0. Solving this linear ordinary diff^erential equation gives 
d{t) = 7ie~''* + 70; where 71 and 70 are two arbitrary constants. Therefore the expressions (|3ip 
and (I36p can be rewritten as 

^ = + 716"''* + 70, /(x) = l, h{x) = b^ mod . 

System (j33P is verified to be true. In addition, rj = 2f Ct+{2Cx+i^)u = 6u— 27i6e~^*. Therefore, 
we have 

S^ = \bx + 7ie~''* + 70, r] = bu- 27166"^*, , , 

/(x) = l, h{x) = b\ ^'''> 

where 6,71,70 are arbitrary constants. Thus, we have a nontrivial regular reduction operator 
Q = dt + i^bx + 7ie-^* + jo)dx + {bu - 2-fibe'^')du. (38) 
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Case 2. 2d: If system (j27|) is satisfied, then the expression (pTj) can be rewritten as (p2|) . 
Substitute it into the first equation of (p3]l to obtain 

b{t){x + q)[h + {x + q)h^] = [b"{t) + 2bit)b'{t)]. 

Substitute the fourth equation of (p7|) into it to get b(t){x + q)[h + {x + q)hx] = 0. It fohows 
that h{x) = r/{x + q), where r is a nonzero constant. Substitute it and (I32p into the second 
equation of 1^ to obtain 2b'^{t) + b"'{t)-rb'{t) = 0. Prom the fourth equation of ([27]), we find 
b"'(t) = 4b^(t)b'{t)-2b'^{t). Substitute it into the preceding equation to get b'{t)[Ab'^(t)-r] = 0, 
which leads to b'{t) = 0. Then a{t) , b{t) , c{t) , d{t) satisfy system (pl|) . and the solution is in- 
cluded in the case 2.2a. 



Case 2.3: When m = 2, system (fT9|l implies 



{7Cx + ^f). = 0, 

2{L + 2ef )6 + + + 4^*. + ^tt = 0, 

2CtCtx + Cut = 0, (39) 
2/(2^0. + C^mx + Cf)t + f{2C. + i^)u - ^hfit - 2/6(2^.. + ef). - 2(/6).x = 0, 

e/ix + 2/iex + (2ex+ef = 0. 

From the first and the last equation of system (j39p . we can get 



f(x) = \^\-'/'exp{a{t)l^), h{x) 
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where a{t) is an arbitrary function, g is a constant. Substituting these expressions into the rest 
equations of system ([5^ . we can see that ^(t,^;) and a{t) satisfy the overdetermined system of 
partial differential equations 



266a: + ittt = 0, 

6t - 36x - ^itix + 4eete + 4a6 + 4at = 0, 

2e{\ix + a)it + 2e2[(i6 + a?]t - eUxx - itiil - + q) 

-4[(at - Utx)ia - |6)e - U%xx] = 0. 

In addition, we have 



(40) 



r? = 2/6 + (26+Cy)^^ = 26ier'/'exp(a(t) J ^) + [i^. + a(t)]n. 

Thus, we have a nontrivial regular reduction operator 

Q = + e(t,x)9, + {26|er'/"exp(a(t) / y) + [^6 + a(t)]n}5„, (41) 

where ^(t,x) and a(t) satisfy the overdetermined system of partial differential equations (j40p . 
In particular, if 6 = 0, from system (I39p we can obtain 

^ = ^(x), ?7 = i(^^+a)u, 

dx, , e-^^^xx+q (42) 



/(x) = |e|-^/^exp(^|^), M^) 
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where a, are arbitrary constants. Thus, we have a nontrivial regular reduction operator 

Q = dt + ax)dx + i\Cx + a)udu, (43) 
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which is equivalent to operator (jlSp with n = 1. Therefore, this special case can be included in 
case 1.2 and we can impose an additional constraint 7^ on the regular reduction operator 

Case 2.4: When m = 3, the third equation of (fTOl) implies 

'^CtCtx + im = 0, 

2/(2ex + i^m. + i^)t + fm. + - 2/6(2^. + c^)x - 2{f^tu = 0, 

6/i/et + (2ex+c^)xx = o, ^^^^ 

^hlf- - - 2h{2ix + ^if) = 0, 
the fourth equation of which can be rewritten as 
1 / fx hx \ 

Since / and h are independent of t, integrate both sides of the preceding equation to obtain 
^ = r(t)|//i|-V4, where r{t) is a function of t. Substituting it into the first equation of (fli]) 
yields the fact that r"'{t) = q{x)r''^{t), where q{x) = -2(|//i|~V^)^. It follows that r'{t) = 
or q'{x) = 0, r"'{t) = qr'^{t). 

Case 2.4a: If r'{t) = 0, then ^ = r|//i|^^/^ (r = const), = 0, = 0, the second equation 
of (jm is an identity, so is the second equation of (fT9|) . The third equation of (jli]) reduces to 
(2^-j; = 0. Combine it with the first equation of ()19p . and use a progress similar to the 

case m = (i.e. 2.1a), we get the expression of ^ and f{x) as stated in ([25]) . where a, 6, c, d 

are constants. From ^ = r|//i|~^/^, we see that h{x) = iby(-)^, where > 0. In addition, 
r] = 2fit + (2^^ + ^^)n = (ax + h)u. Thus, we have 

= 2ax^ + 46x - cx + d, 
rj = {ax + 6)n, 



2ax^ + 46x — cx + d 

where a, 6, c, d, r are arbitrary constants. Thus, the corresponding regular reduction operator 
has the form 

Q = dt + (2ax^ + 46x — cx + d)dx + (ax + b)udu, 

which is equivalent to Lie symmetry operator. 

Case 2.4b: liq'{x) = 0,r"'{t) = qr'^{t), then (|//i|"^/'')x = -^q- Integration both sides of it 
gives = —\qx + s, where s is a constant. Therefore ^ = r{t)\fh\~^^^ = r{t){—^qx + s). 

From the first equation of ([IHD we see that (2^^; + C^)^ = jCxx + 3(7^^ + 4Cy)x = \Cxx- 
Substituting the last two expressions into the third equation of (j44p . yields r'(t) = or ^ = 0, 
which are the cases have already been discussed. 



Case 2.5: When m 7^ 0, 1, 2, 3, the third equation of (|T9|) implies 
^hff- - + (1 - m)h{2ix + i^) = 0, 

6 = 0, (46) 

{2^x C~J')xx — 0. 
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Notice that the second equation of (06]) indicates that ^ is independent of t, therefore the second 
equation of (jl9p is satisfied automaticahy. Similar to the case of m = 0, from the third equation 
of ([Ml) and the first equation of we get the expression of ^ and f{x) as stated in ([28j) . where 
a, b, c, d are constants. Substituting the expression of ^ and / into the first equation of (06]) we 
obtain 

, / f -7ax - 7b + 2c + (1 - m)(ax + b) , \ 

h(x) = exp { / 7i r, ; ax). 

\ J 2ax^ + 4bx - cx + d J 

In addition, t] = + (2^a; + = (ax + b)u. Therefore, we have 

^ = 2ax^ + 4:bx — cx + d, 
7] = {ax + b)u, 

/W-e.p(/- „J°-^-^''t^' d.). (47) 



2ax2 + Abx — cx + d 

, f" f -7ax -7b + 2c+ il - m){ax + b) , \ 

= exp / — 5 ; dx . 

^ ' \ J 2ax^ +4:bx- cx + d J 

where a, b, c, d are arbitrary constants. Thus, the corresponding regular reduction operator has 
the form 

Q = dt + (2ax^ + ibx — cx + d)dx + {ax + b)udu, 
which is equivalent to Lie symmetry operator. 



From the above discussion, we can arrival at the following two theorems. 

Theorem 8. A complete list of -inequivalent equations ([3]) having nontrivial regular reduction 
operator is exhausted by ones given in tablelM 



Table 2. Results of regular reduction operator classification of class (|3]) 



N 


n 


m 


/(^) 


h{x) 


Regular reduction operator Q 


1 


/-I 


71+1 


^ pfSrl exp (r / 1 dx) 






4(n + l)C2 


2 


1 


1 


1 


dt + i^bx + 716"''* + 7o)& + {bu - 2-/ibe-'")du 


3 


1 


2 


|C|-5exp(«(t)/f ) 


85^ 


dt+at,x)d. + {26[er 5 exp(a(t) / f ) 
+ + a{t)]u}du 



Here r,p, 6, 71,70 are arbitrary constants, ^{x) in case 2.1 is an arbitrary functions of the variables x, ^{t,x) and 
a{t) in case 2.3 satisfy the overdetermined system of partial differential equations (140 |l and ft 7^ 0. 

Theorem 9. Any reduction operator of an equations from class ([3]) having the form ()12p with 
Ct = 0, ^xxx = is equivalent to a Lie symmetry operator of this equation. 



4 Exact solutions 

In this section, we construct nonclassical reduction and exact solutions for the classification 
models in table [2] by using the corresponding regular reduction operator. Lie reduction and 
exact solutions of equation from class ([3]) have been investigated in reference |15j . We choose 
case [J in table [2] as an example to implement the reduction, the other cases can be considered 
in a similar way. 

For the first case in table [2l the corresponding equation is 

[I exp (r 1 i dx)]uu - {u^u,), - ^^^(^^ffj^^"^' = 0, (48) 
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which admit the regular reduction operator 

Q = dt+ ^{x)d^ + {- + ,r^)ndu. 

n 2n + 2 

An ansatze constructed by this operator has the form 

1 r f 1 f 1 

x) = (/5(a;)| ^1 2n+2 exp (— / -dx), where u = t — / - dx. 

n J i ' J ^ 

Substituting this ansatze into equation (j^H]) leads to the reduced ODE 

[(4r2 -p)n2 + 4(2n + l)r'^](p''+'^ (uj) + 4n(n + l)[n(f"{uj) - 2(n + l)r(p'{uj)]ip''{u}) 
+4n3(n+l)99'2(tj)(^«-i(tj)-4n2(n+l)v9"H =0. ^ ^ 

Because there are higher nonlinear terms, we were not able to completely solve the above equa- 
tion. Thus, we try to solve this equation under different additional constraints imposed on p 
and r. 

We first rewrite equation (jl9|) as 



4n2(n + l)[>f'{uj)(p''{u;)]' - 4n^{n + l)ip"(u;) - 8n(n + l)'^rip'{uj)(p''{u;) 

+ [4(n + l)2r2 - pn2](p"+i(a;) = 0. ^ ' 

If we take p = 4(1 + ^)'^r'^, then the general solution of ()50p can be written in the implicit 
form 

Up to similarity of solutions of equation ([3]), the constant C2 is inessential and can be set to 
equal zero by a translation of cv, which is always induced by a translation of t. 

If we further set n = 1, the general solution (|5ip can be rewritten in the implicit form 

ln(2V_^^) ^arctanh(^) 

w + C2 = 0. (52) 

4r 2^/cir 



Thus we obtain the following solution 

u{t,x) = ip{u})\£,\^ exj)(r [ ^ dx) 



cjj = t — I - dx 



for the equation 



[I ^1 I exp J ^ dx)]utt - {uu^)^ - — — — —u^ = 0, 

where (p satisfy the equation (j52p . ^ is an arbitrary function and r is a non-zero constant. 
If we further set r = 0, the general solution ()5ip can be rewritten in the implicit form 

^ (^"+i(a;) -V9(w) = cia; + C2. (53) 



n + 1 

Thus we obtain the following solution 

u{t,x) = ip{u})\(\^, uj = t-J^ 



dx 
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for the equation 

where (p satisfy the equation ()53p . ^ is an arbitrary function. In particular, for n = 1 from 
equation (f53]l we have 



ifioj) = 1 ± ^/l + 2{clU} + C2). 
Thus we obtain an exphcit solution 

u{t,x) = [1± Vl + 2(cit^ + C2)]|C|3, u} = t-j^dx 

for the equation 

I 41 iUtt-{uUx)x 

If we take different functions for ^, then we can obtain a series of solutions for the corresponding 
equations. In order to avoid tediousness, we do not make a further discussion here. 



5 Conclusion and Discussion 

In this paper we have given a detailed investigation of the reduction operators of the variable 
coefficient nonlinear wave equations ([1]) (equivalently to ([3])) by using the singular reduction 
operator and the equivalence transformation theory. A classification of regular reduction oper- 
ators is performed with respect to generalized extended equivalence groups. The main results 
on classification for the equation ^ are collected in table [2] where we list three inequivalent 
cases with the corresponding regular reduction operators. Nonclassical symmetry reduction of 
a class nonlinear wave model (I48p which are singled out the classification models are also per- 
formed. This enabled to obtain some non-Lie exact solutions which are invariant under certain 
conditional symmetries for the corresponding model. 

The present paper is a preliminary nonclassical symmetry analysis of the class of hyperbolic 
type nonlinear partial differential equations ([T|). Therefore, further investigations of different 
properties such as nonclassical potential symmetries, and nonclassical potential exact solutions 
as well as physical application of this class of equations would be extremely interesting. These 
results will be reported in subsequent publications. 
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